Kink interactions in SU{N) x Z2 
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There are A'^ — 1 classes of kink solutions in SU(N) x Z2. We show how interactions between 
various kinks depend on the classes of individual kinks as well as on their orientations with respect 
to each other in the internal space. In particular, we find that the attractive or repulsive nature 
of the interaction depends on the trace of the product of charges of the two kinks. We calculate 
the interaction potential for all combinations of kinks and antikinks in SU{5) x Z2 and study their 
collisions. The outcome of kink-antikink collisions, as expected from previous studies, is sensitive 
to their initial relative velocity. We find that heavier kinks tend to break up into lighter ones, 
while interactions between the lightest kinks and antikinks in this model can be repulsive as well as 
attractive. 



I. INTRODUCTION 

Topological defects are observed in condensed matter 
systems and may have been formed during phase transi- 
tions^ at early stages in the history of the universe [Q. 
If observed, they virould provide invaluable information 
about the universe when it was a tiny fraction of a sec- 
ond old. If not observed, topological defects still play an 
important role by placing constraints on particle physics 
models and cosmology. The formation and scaling of a 
network of defects strongly depends on how they interact 
among themselves. This, in turn, will affect the type and 
the strength of restrictions that observations (or the lack 
thereof) can impose on the underlying model. Another 
context, in which interactions between defects are impor- 
tant, is the possible connection between elementary par- 
ticles and solitonic solutions of classical field equations 

i- 

In addition to the more general reasons given above, 
understanding how the SU{N) x Z2 domain walls interact 
could be important in the Hght of the correspondence, 
found by Vachaspati [D, between the spectrum of SU{5) 
monopoles and the spectrum of one family of fermions in 
the Standard Model. Interactions between SU{N) x Z2 
kinks may also be relevant to the solution of the monopole 
over-abundance problem based on sweeping monopoles 
with domain walls as proposed in |^. 

Previous work on interactions between kinks has 
mainly concentrated on the sine-Gordon and the (j>'^ mod- 
els in (1+1) dimensions ||, |, |, |, |, 0, 0, |l|]. Even a 
relatively simple system of a kink interacting with an 
antikink can have rather non-trivial dynamics, which is 
one of the reasons why so many researchers have worked 
on this problem in the past. The force between kinks and 
antikinks of the 0* model is always attractive. The out- 
come of their collision can be one of the three types: they 
can annihilate, they can scatter off each other or they can 
form an intermediate bound state before ultimately sep- 



arating or annihilating. The general tendency is that at 
low collision velocities kinks tend to annihilate and at 
higher velocities they scatter. However, the dependence 
of the outcome on the incident velocity is rather non- 
linear, as was found in ^ |^, |^, ^, ||, |l^ and investigated 
in detail in jl^. Namely, it was found that, over a rela- 
tively small range of initial velocities, intervals of initial 
velocity for which kink and antikink capture each other 
alternate with regions for which the interaction concludes 
with escape to infinite separations. In this alterna- 
tion phenomenon was attributed to a nonlinear resonance 
between the orbital frequency of the bound kink-antikink 
pair and the frequency of characteristic small oscillations 
of the field localized at the moving kink and antikink 
centers. 

In this work, when discussing interactions between 
kinks, we will aim to concentrate on issues that are 
unique to SU{N) x Z2 and will refer to earlier work 
when a problem can be reduced to that of kinks in the (j)'^ 
model. In particular, we will show that in SU{N) x Z2 
kinks can repel as well as attract. 

This paper is organized as follows. In Sec. || we give a 
brief overview of kink solutions in SU{N) x Z2. In Sec. 
Ill we develop a framework in which SU{N) x Z2 kink 
interactions can be discussed and show a simple way of 
determining whether a given pair of kinks will attract or 
repel. In Sections |l^ and ^ we study the kink-antikink 
interactions in SU{5) x Z2. Results are summarized in 
Sec. 



II. KINKS IN SU{N) X Z2 

Consider a (l+l)-dimensional model of a scalar field 
<& transforming in the adjoint representation of SU{N), 
with N taken to be odd and with the additional Z2 sym- 
metry that takes $ to — $. The Lagrangian is 



^ Here, the term "phase transition" includes continuous transitions 
called crossovers. 



L = Tr(a^$)2 - !/($) 



(1) 
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where V{^) is such that $ has an expectation value 
that can be chosen to be 



nln+i 
N{N^-1) V 





-(n + l)l„ 



(2) 



where Ip is the p x p identity matrix, N = 2n + 1 and 
77 is an energy scale determined by the minima of the 
potential V. Such an expectation value spontaneously 
breaks the symmetry down to: 

H = [SU{n + 1) X SU{n) x C/(l)]/[Z„+i x Z„] . (3) 

Various types of kink solutions in this model are de- 
fined by choices of the boundary conditions at x — —00 
and X — +00, where x is the space coordinate. It was 
proved in that for a kink solution to exist one must 
necessarily have [$(—00), $(+00)] = 0. This allows one 
to Hst all the possible boundary conditions (up to gauge 
rotations) that can lead to kink solutions. We can fix 
<f>L = $(— cx)) = $y given in Eq. (Eh. Then we can have 



$(+00) 



eTV\ 



N{N^ - 1) 



diag(nl„+i_,, -(n + 1)1,, nlg,-{n + l)l„-q) , (4) 

where we have introduced a parameter ex — ±1 and an- 
other g = 0, n. The label et is -1-1 when the boundary 
conditions are topologically trivial and is —1 when they 
are topologically non-trivial, q tells us how many diago- 
nal entries of $l have been permuted in ^fi. The case 
g = is when — ct^l- As was suggested in |Q, 
the lowest energy stable topological {ct = — 1) kink solu- 
tion corresponds to = n. Topological q — n kinks were 
studied in detail in 

The most general ansatz for the kink solution was 
found in ||l^ and can be written as: 



$fc = F+{x)M+ + F_(a;)M_ + gix)M 



where 



M+ = 



M = 



(5) 



(6) 



and, for q ^ and q ^ n, 



M = ndiag{q{n - q)ln+i-q,-{n - q){n + 1 - q)lq, 

-{n - q){n + 1 - q)lq,q{n + 1 - q)ln-q) (7) 

with 

^i^^^[2q{n~q){n + l-q){2n{n+l-q)-q}]-^/^ . (8) 

For g = or for q = n, the matrix M is zero. The 
boundary conditions for F± and g{x) are: 

F^{±oo) = ±l, F+(±oo) = l, g{±oo) = 0. (9) 



One can define the charge of the kinks as 

g= i($fc(-h00)-$fe(-C50)), 

which corresponds to a current 



1 



(10) 



(11) 



where = 0,1, and e'^" is the antisymmetric tensor. 
The definition (|lO|) of Q can be used for non-topological 
kinks as well as topological ones. 



III. KINK INTERACTIONS IN SU{N) x Z2 

Consider two kinks, if'-^' and K'^^\ separated by a dis- 
tance which is larger than their core sizes. The classes to 
which the two kinks belong, as well as the global topol- 
ogy of the two-kink configuration, are determined by the 
choices of the three vacua: 

• $_ at X = —00, 

• $0 in between the two kinks, 

• $+ at X = +00 . 

Let indices {ei^\q^^'') describe the kink between <&_ 
and $0, where g*-^-' denotes the kink class as defined in 
Sec. 0, and e^^^ is —1 for topological and -fl for non- 
topological kinks. Similarly, let (e^\g*^^^) denote the 
kink bounded by $0 and Boundaries $_ and $_)- of 
the two-kink system can also be described using indices 
(e^^ , q^^^ ) defined in a similar way . Topology requires 



that e 



(3) 



Thus, the two-kink system can be 



described by (q^^\ q^'^\ q^^\ ei^\ eip^ ) . This notation is 
invariant under global gauge rotations, since it contains 
only information about how many diagonal entries were 
permuted in each of the vacua $0 and $+ with re- 
spect to each other. 

For given values of g'-^' and g'^^^ , not all values of g'^' 
will generally be allowed. To determine the selection pro- 
cedure, let us start with $0, in which g^^-* diagonal entries 
of <&_ were permuted. We want to know how many diag- 
onal entries of $_ can be permuted in <i>+ (i.e. the value 
of 9^^^) given that g^^^ diagonal entries of $0 were per- 
muted in There are n+1 entries with absolute values 
equal to n in $_ (see Eq. (||)), which we can denote by 
^'s, and n entries with absolute values equal to (n + l), 
which we will denote by B^s. We will refer to those ^'s 
and S's of $0, that were permuted to form $0 out of 
as "changed", and the ones that were left untouched as 
"unchanged". Then, when permuting ^'s and S's of $0 
to form <f>+ (q'--^^ permutations ), one has the following 
options: 
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1. Permute a changed A with a changed B. This op- 
eration will decrease the value of g^^^ by 1, since 
it will restore the original order of the given pair 

(2) 

of A and B in Let q\ denote the number of 
possible ways it can be done. Since there are only 
q^^^ changed A's, gf-* < min(g(^), g^^)). Also, if 
there is a deficit of unchanged ^'s, one is forced to 
permute at least (g*^^-' — n + q^^^ — 1) changed A's, 
which means that q^'^ > max(0, g*^^' — ri + g'^' — 1). 

2. Permute a changed A with an unchanged B. This 
operation does not affect the value of g^'^^. The 
number of possible ways in which this can be 

(2) 

done, denoted by 53 , is Hmited by the num- 
ber of available changed A^s and unchanged S's: 

(2 
I2 



-^'^ <min((z(i)-(?f\n-(7(i)). 



3. Permute an unchanged A with a changed B. This 
operation also does not change the value of q'^^^ . It 
can be done in gg ways, limited by the number of 



available unchanged ^'s and changed S's: q^ 



(2) 



< 



min(ri, — q'^^") + 1,9 



(2)n 



4. Permute an unchanged A with an unchanged B. 

This operation will increase the value of q'^^^ by 1, 

(2) 

and can be performed in q\ ways, limited by the 
number of available unchanged A^s and unchanged 
B's: qf^ < mm{n - - q^ + l,n~ 



In summary, independent of and , given q^^^ and 
q^'^\ the set of possible values of g'^^ can be found using 



(3) 



(1) 



(2) 



14 



(2) 

where we have defined integers ql {i — 1,2,3,4) that 
can take on all non-negative values allowed by the fol- 
lowing selection rules: 

max(0,9(^^ - n + g'^' - 1) < < min(g(i) , g'^)) ^ 
0<g^^' <min(q(i) , 
< gf^ < min(n-- + - gf^) , 

< q^i^ < min(n - q^^'^ - qf^ + 1, n - q^^^ - q!^^) .(13) 

Next we would like to determine whether a given pair 
of kinks will attract or repel. We can use a procedure 
which is analogous to one used by Manton ||lj| in the 
case of the 0^. At large separations, the two-kink ansatz 
can be written as: 

$(a;) = {x + a) + (a; - a) - $0 , (14) 



(2) 



(12) 



where ^^i}\x) is the first kink solution with $^^^(—00) = 

(2.) ('2I 

$_ , $^ (a;) is the second kink solution with $j, (-I-00) = 
$+, $0 = *i^''(+oo) = $^^•'(-00) is the vacuum in the 



(1) 



region separating the kinks and a > is the distance 
from the origin to the centers of the kinks. Using Eqns. 
(^ and (|To|), we can write: 

In ||l^ it was shown that functions F^'^^ and g^^'^^ can be 
treated as approximately constant for a relatively wide 
range of parameters of a general quartic potential. A 
simplified (and, therefore, only approximate) version of 
the two-kink ansatz is given by: 



$(x) « Af|^) + l^l^^Qd) 



where 



^ + $0 



(16) 



(17) 



and 



w tanh[cr(a; + a)] , w tanh[cr(x - a)] , (18) 

where a ^ and <7~^ is the "width" of the wall. 

The field energy-momentum tensor can be derived us- 
ing the action principle and is given by 

rj^f.,^ ^ 2Tr[(9^$a''$]-r;^''Tr[9'"$a„$]-r;^'^y($) . (19) 
Therefore, the momentum density in (1+1) dimensions 



IS 



7^ = = -2Tr[$$'] 



(20) 



where $ = 9f$ and <&' = dx^- One can define the mo- 
mentum of a field configuration on the interval Xi < x < 

X2 as 



P 



2Tr[$$']da 



(21) 



To calculate the force between the kinks, let us consider 
the initial rate of change of momentum of a (intially) 
static field configuration given by the two-kink ansatz: 



dP 



2(Tr[$$'] + Tr[$$'])dx 



(22) 



One can use field equations of motion and integrate to 
obtain 



dP 
'dt 



Tr[$^] - Tr[($') ] + ^(*) 



(23) 



Let us choose Xi ^ —a and — a <C 0:2 <C a (e. 5. 0:2 = 0). 
That is, we want to estimate the force on the first kink 
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(the one aX x = —a), due to the second kink. Let us 
define 

/ = tanh[cr(x + a)] and x = tanh[o'(a; — a)] + 1. (24) 
The two-kink ansatz given by (|l5|) can then be re- written 



H^)^^^ + lQ^'^f + lQ^'\~i + x). (25) 

Initially, $ = 0. Also, within the range xi < x < X2, 
x{x) ^ 1 and we can perform an expansion in x- To the 
leading order in x we find: 



dP 



^Tr[(QW)2](/')2 - ^Tr[Q(i)g(2)]/V 



dV - 



[^] x=0 + X! [ 5$a ] x=0 ^ 



where functions x"' are defined by 
and coefficients [Q^^-*]" are defined by 



(26) 



(27) 



(28) 



where are the SU {N) generators normalized so that 
Tt\T°-T^] ~ 5ab/'2.. Using the equations of motion gives 



dP_ 



Tr[Q«Q(2)]/Y + ^[$'^"]x=o 



X2 
Xi 



= ^Tr[QWQ(^)][-/Y + /"xl" 



(29) 



It can be shown that 



f'x' + f'x 



< for all xi and 



X2 satisfying the constraints. Thus, the sign of dP/dt, 
and the attractive or repulsive nature of the force, is de- 
termined by the sign of Tt[Q'-^'^Q'^^^. 

Next we would like to express Tt[Q'-^'> Q'^^'>] in terms 
of parameters {q^^\ q^'^\ q''^\ e^^) which define the 
two-kink configuration. We write 



Tr[Q(i)Q(2)j ^ _(Tr[$_$o]+Tr[$o*+] 
- Tr[$_$+]-Tr[$o*o]). 



(30) 



Applying definitions (||) and (Q) to boundary conditions 
specified by {^l, ^r) we find: 



2eT 



N{N^ - 1) 
2q''^\n + l)n + {n - q){n + Vf] 



2 



1-g 



(2n+l) 
n(n + 1) 



(31) 



Using the above expression and the fact that Tr[$o^o] 
?7^/2 we find: 



Tr[Q«Q(^)] = i(4^) 



(2) _ (1) (2) 



1) 



One can see, for example, that in the case of a topological 
kink interacting with a topological antikink (e^"* = — 1 

(2) 

and = ~1) both, attraction and repulsion, are possi- 
ble depending on the choices of g^^-' , g'-^-' and g'-^-' . This is 
a novel feature, when compared to the classical (f)^ case 
1^, where the force between a kink and an antikink is 
always attractive. An analogous situation is found in the 
case of interactions between global 0(3) monopoles 
There, attraction or repulsion between a monopole and 
an antimonopole is determined not by their topological 
charges but by the relative phase of their field configura- 
tions. 

Our derivation of the fact that the sign of Tr[(3(^^Q(^)] 
determines whether kinks will attract or repel was inde- 
pendent of a particular form of T^($)- We did, however, 
rely on certain approximations in the derivation that may 
not be valid for some extreme choices of V{<^). Our nu- 
merical investigation of kink interactions in SU{5) x Z2 
with a general quartic potential has always yielded an 
agreement between the sign of TtIQ^^^^ Q'^^^] and the over- 
all sign of the interaction potential. 

In the next section we consider the = 5 case and 
study interactions between kinks and antikinks in more 
detail. 



IV. KINK INTERACTIONS IN SU{5) x Z2 

Let us consider the (l+l)-dimensional model of a 
scalar field $ transforming in the adjoint representation 
of SU{5) X Z2. We will take the potential to be 

y($) = -m2Tr[$2] + /i(Tr[$2])2 + ATr[$'^] + Vo , (33) 

with parameters such that the vacuum expectation value 
(VEV) of <& breaks the symmetry spontaneously to 
[SU{3) X SU{2) X C/(l)]/[^3 X Z2]. This happens in the 
parameter range 



The VEV can be chosen to be <i>y 
7? diag(2, 2, 2, -3, -3)/(2VT5), where 



and 



(34) 



(35) 



(36) 
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Tr[Q(^)Q(2)] for 


= 


g(^) = 1 


gf^*) = 2 


qW = 0, = 


-2 (A) 


X 


X 


qW = 0, q'-^> = 1 


X 


-1 (A) 


X 


= 0, = 2 


X 


X 


-k (A) 


q(l) ^ 1^ g(2) ^ 1 




-! (A) 


-1 (A) 


g(i) = 1^ g(2) ^ 2 


X 


(A) 


+ T2 (R) 


q(i) = 2, gt^) = 2 


-1 (A) 


+ ^ (R) 


X 



Table I: Tr[g(i'Q(2'] evaluated for different choices of the 
topological kink - topological antikink configuration in SU{5). 
(A) denotes attractions, (R) - repulsion and x - means that 
the arrangement is impossible. 



The constant Vb — m'^Tf'/A in Eq. (33) ensures that 
V{^v) = 0. 

We will study interactions between topological kinks 



and antikinks 



-1, e 



(2) 



— 1) with all possi- 
ble choices for indices g^^', g*-"^^ and q'^'^\ as defined in 
Section |l|. We will not consider interactions between 
non-topological kinks nor interactions of non-topological 
kinks with topological ones. Our initial configuration will 
be that of two well-separated solitons moving towards 
each other. From here on, we will label such configura- 



tions with three indices. 



assumed that e 



(1) 



,(2) 



(^(1), ^(2) 
= -1. 



'''^-'), and it will be 



In this model, only the q = 2 topological kink solution 
is stable. It also has the smallest energy of all kinks. 
However, we will not restrict the analysis to g = 2 kinks, 
as interactions between all types of kinks antikinks could 
be potentially interesting. 

We start by evaluating Tr[Q^^^(5*-^''] for all possible 
choices of the configuration of a kink and an antikink 
using Eq. (^. As we have shown in Section III, this 



whould tells us whether the two kinks will attract or re- 
pel. The results are given in Table | and show that there 
are two possible combinations in which there is a repul- 
sive force between the solitons. 

We would like to evaluate the interaction energy be- 
tween different classes of kinks and antikinks in SU(h) x 
Z2. The cases when analytical kink solutions are known 
are the q — Q kink and, for the special value h/\ = 
-3/20, the g = 2 kink The g = kink solution 

is the same as in the (j)^ model, since in this case the only 
non-zero component of $ is the one along $y, that is 
<I> = (f){x)^v /"H, and the potential takes the same form as 
in the (f)^ theory: 



(37) 



Therefore, the interaction potential between q = kinks 
and antikinks will be identical to the one found in the 0'* 
model 0. 

Let us next consider q = 2 kinks and antikinks. Con- 
sider a kink at a; = — a and an antikink at a; = a, with 



a > 0, each moving with a velocity v directed towards 
the origin. For values of a larger than the core sizes of 
two kinks the following ansatz is valid: 



KKix)^^Kix + a) + ^Kix-a)-^o, (38) 



where $0 = ^k{+oo) = ^j^{—oo) is the field in between 
the two kinks. For h = — 3A/20 and q — 2 the kink 
solution is known flla, M: 



*,=2(-^) = 



[1 — tanh(o'a;)] ^ ^ I"*" ^ tanh((Tx)] 



(39) 



where a — m/V2. The ansatz for the kink and antikink 
can then be written as 



2 
2 



[1 + Fk] 
2 

[1 + Fk] 



(40) 



where $_ = $^^(-00), $0 = *Ki?(0), = 
$Kif(+oo), 



Fk — tanh[cr7(x + a)], 
Fg- = teL'Dh[(j"f(x — a)], 



(41) 
(42) 



and 7 = l/Vl — v'^ is the Lorentz factor. 

The total energy of the ansatz ( |4Q| ) is a sum of the gra- 
dient and the potential energies obtained by integrating 
corresponding energy densities along the space direction 
x: 



E = G + P, 



where 



G 



+ 00 



dx Trid^^KR) 



(43) 



(44) 



r+°° 

P= / dx{-m'TTi^KKf + h{TTi^KKfy 

J —00 



2^2 



■} (45) 



At first let us consider a q = 2 kink and a. q — 2 
antikink such that $i^(-oo) = ^izi+oo) (the g'^' = 2, 



(2) 



2, q 



(3) 



case in Table mj. A possible set of 



boundary conditions corresponding to this case is: 



$0 = 



2V15 



2V15 



2V15 



diag(2,2,2,-3, -3), 
diag(3,3,-2,-2,-2), 
diag(2,2,2,-3,-3). (46) 
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Using these boundary conditions in Eq.(pO|) and substi- 
tuting latter into Eqns. and (|45| ) gives 



r + oo 3 

G = y dX-^idxFK-dxF^ 



K)^ : 



- 4(^^K - + A{Fk - F^f] , (47) 

where X — ^mx/\/2 and we have taken h/X ^ —3/20. 
Evaluating integrals in Eq.(p7|) and using R = V2j ma 
yields 



G = 



A ^3 



sinh R — R cosh _Rn 



sinh-^ R 



-1 



1 



sinh^^i?^ 4 



R) 



-«(Zi?-l) + e^(ii?-l)}] 



(48) 



As expected, at i? = oo the total energy is equal to the 
sum of the two kink masses (2x4V2m3/(3A))||l|| divided 
by the Lorentz factor. Subtracting E{oo) from E{R) = 
G{R) + P{R) leaves only the interaction part of the total 
energy: 



UoiR) = 



A sinh^ R 



sinh R - 



,1 



R cosh R 
1, 



(49) 



The dependence of Uo{R) on i? for 7 = 1 is shown as 
a soHd Hne in Fig. 0. It clearly indicates an attraction 
between the kink and the antikink. 

The validity of the ansatz ( ^8|) cannot be justified for 
small values of R. To test the analytical result, we have 
evaluated the interaction energy numerically by explicitly 
integrating the full set of equations of motion and evalu- 
ating the spatial integral over the gradient and potential 
energy densities at each time step. The separation R be- 
tween the kinks was defined as the distance between the 
maxima of their potential energy densities. The shape 
of the found interaction potential was similar to the one 
given by Eq. (||), not only for h/X = -3/20, but for all 
considered values of the parameter: —7/30 < h/X < 100. 
The interaction potential for walls with different values 
of 7 > 1 was also qualitatively the same. 

Next, let us consider the interaction of a g = 2 kink and 
a q = 2 antikink but with $/<-(— 00) 7^ ^j^{+oo). This 
would correspond to the {q^^^ — 2, g^^^ — 2, q^^^ = 1) 
case in Table |. Keeping and $0 the same as in Eq. 
(H) we can choose to be 



V 



2V15 



diag(2,2,-3,-3,2) 



(50) 




Figure 1: Interaction potentials Uo, given by Eq. ( |4s| ) (solid 
line), and Ui, given by Eq. dsi] ) (dotted line), for 7 = m = 
A = 1. 



This choice of boundary conditions leads to 

G = J dX^{2dxF^ + d^FKdxFji + 2d,F^) 



P = 



+ 00 



dX- 



3^-[8(F.-^,)^ + 8(^^.-F,)3 

U{Fk - Fj,f - 30(Fk - F^)(l - FkFi,) 
mFKF^iFl + F| - 2) + 10Fa'F^(1 - F^^F^) 



35(1 - FlFl)] 



(51) 



Performing the integration and subtracting E{<X)) gives 
the interaction potential between the two walls: 



Ui{R) 
1 

^8^ 



^/2r"^ 



m 7 



A sinh'' R 



R cosh R — sinh R 



[e-3«(3 + 7i?) _ e-'"(ll + 13i?) 

-Ke^(8-4i?)]) 



(52) 



We find that in this case, the interaction potential is 
purely repulsive. The dependence of Ui(FC) on R for 
7 = 1 is shown as a dashed line in Fig. |l|. The ex- 
act numerical evaluation of the interaction energy of this 
kink-antikink system, using the method outlined above, 
did not show qualitative deviations from Ui{R) for all 
considered values oi h/X and 7. 



The two configurations: {q^^ 
and (g(i) = 2, = 2, ^ 
studied (^''-equivalent case (g*^^-' 



) = 2. 



!, (7(3) = 0) 
1), as well as the well- 
= 0, = 0, = 0), 
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are the only ones for which analytical kink solutions are 
known. Other configurations from Table | were treated 
only numerically. For all considered choices of parame- 
ters and velocities, we did not see any deviation from the 
redictions for the attraction or repulsion given in Table 



V. SU{5) X Z2 KINK-ANTIKINK COLLISIONS 

In this section we will study kink-antikink collisions. 
Even in the "simple" case of a (j)'^ kink colliding with an 
antikink the variety of possible outcomes is surprisingly 
rich d, |6[ I, I, §, ^ Depending on the incident ve- 
locity, the two kinks can annihilate, or they can bounce 
off each other and never meet again, or they can form 
an intermediate bound state, namely, they can bounce 
off each other several times before ultimately separating 
or annihilating. The dependence on the incident velocity 
is rather non-trivial, as was found in g I, |, I |, HI 
and investigated in detail in jl^. Namely, it was found 
that, over a relatively small range of initial velocities, 
intervals of initial velocity for which kink and antikink 
capture each other alternate with regions for which the 
interaction concludes with escape to infinite separations. 
In this alternation phenomenon was attributed to a 
nonlinear resonance between the orbital frequency of the 
bound kink-antikink pair and the frequency of character- 
istic small oscillations of the field locaHzed at the moving 
kink and antikink centers. 

We will not attempt a study of exact dependence of 
outcomes of kink-antikink collisions on initial velocities. 
The reason is that in SU{5) x Z2 there are too many 
possible combinations and there is an additional param- 
eter, h/X, which can possibly affect the stability of kink 
solutions and the outcome of kink-antikink collisions. In- 
stead, we will describe the outcomes for each of the com- 
binations Hsted in Table | and illustrate the most inter- 
esting ones^. 

In order to study the collisions, we need to integrate the 
field equations of motion forward in time. Without loss 
of generality, we can choose the initial kink-antikink con- 
figuration to be diagonal Since equations of motion 
preserve the diagonal form, one only needs to consider 
the evolution of four functions: a{x,t), b{x, t), c{x, t) and 
d{x, t), defined as: 

$(a;,t) = a{x,t)X3+b{x,t)Xs+c{x,t)T3+d{x,t)Y , (53) 

where A3, As, T3 and Y are the diagonal generators of 
SU(5): 

A3 = idiag(l,-l,0,0,0) , 



^ Several animated kink collisions can be 
^ttp:/ /tlieory.ic.ac.uk/~LEP/su5kinks.htm^ . 
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Figure 2: Kink-antikink collision in the (g^^-* = 0, g'^-* 



0, 



_ g-j gj^gg^ ill ^iig parameter range when q = kink is un- 
stable. The solid line shows the initial energy density profile 
and the dotted line - the final. The right moving q = kink 
collapsed into & q — 2 kink going to the left and the remainder 
q = 1 non-topological (et = 1) kink moving to the right. The 
originally left moving q = antikink splits into a right mov- 
ing q = 2 antikink and a left moving q = 1 non-topological 
kink. The arrows show directions and approximate relative 
magnitudes of kink velocities. 



As = 

T3 = 
Y = 



2%/3 



diag(l, 1,-2, 0,0), 



-diag(0, 0,0, 1,-1) 



1 



2V15 



diag(2,2,2,-3,-3) 



(54) 



viewed at 



As mentioned in Section III , the q — wall is identical 
to the kink of the simplest t^^model. The q = solution, 
however, is known to be locally unstable against pertur- 
bations along diagonal components of $ in the parameter 
range h/X > -3/20 ||||. The outcome of the g = kink- 
antikink coUision will, therefore, depend on h/X as well as 
on Vinitiai- In Fig. ^ and Fig. ^ we illustrate the coUision 
for Vinitiai = 0.05 and h/X — 0. Just before the collision, 
both the kink and the antikink collapse. More detailed 
analysis of functions a, 6, c and d in Fig. || reveals that 
the q = kink, initially moving to the right, has collapsed 
into a g = 2 kink, moving to the left, and the remainder 
q — 1 non-topological (et = 1) kink, moving to the right. 
The originally left moving q = antikink has split into 
a right moving q = 2 antikink and a left moving q = 1 
non-topological kink. The q = 2 kinks will separate to 
infinities, while the fate of q — 1 non-topological kinks 
needs more explanation. As was found in iQ, the q = 1 
non-topological kink is unstable against a collapse into a 
pair of g = 2 kinks. Therefore, depending on the values of 
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-100 -50 



50 100 



Figure 3: Functions a{x,t) (dotted line), b{x,t) (short dash 
line), c(x,t) (long dash line) and d(x,t) (solid line), defined 
m Eq. (M) at the same final snapshot as in Fig. ^ 



Vinitial and h/X, the non-topological q — 1 kinks (in the 
center of Fig. ||) can either immediately decay into radi- 
ation, or they can spUt into pairs of q = 2 kinks moving 
away from each other. We found that the latter was the 
case for the choice of parameters corresponding to Fig. 
I and Fig. |. For -3/20 < h/X < -3/70, when q = 
kinks are locally stable, outcomes of their coUisions are 
the same as in the case of kinks and antikinks in cj)'^ model, 
the case studied extensively in §, |, 0, |, ||, |l|| . 

Fig. ^ illustrates a coUision of a g = 1 kink (initially 
on the left) with & q = antikink (initially on the right). 
The parameters are Vinuiai — 0.2, rj = 1, X — 0.5 and 
h/X = —0.1. The q = kink is unstable for these pa- 
rameters. The mass of the q — kink is 2\^m^/X' ^ 
and, for our choice of parameters, is equal to 0.243. The 
q = 1 kink mass can only be found numerically and is 
0.150. For comparison, the mass of the q = 2 kink with 
these parameters would be 0.033 - almost one fifth the 
mass of the q = 1 kink. We find that, during the colli- 
sion, the q — kink collapses into a q = 2 kink traveling 
to the right and a. q = 1 non-topological kink traveling 
to the left. At the same time, the original q = 1 kink 
collapses into three q — 2 kinks, with outer kinks moving 
away from each other. The final configuration is that of 
four q = 2 kinks arranged so that the two left most kinks 
and the two right most kinks form (g*^^^ = 2, q^^'' = 2, 
= 1) combinations, while the two inner kinks form a 
(q(i) = 2, = 2, = q) combination. The rest of 
the original energy is radiated away. 

In Fig. ||the initial configuration is (g*-^-* = 0, g*^^^ — 2, 
=2). In this case, the q = antikink (initially 
on the right) splits into a q = 2 topological antikink, 
which starts interacting with the q = 2 kink, and a q = 1 



-100 -50 50 100 



Figure 4: Kink-antikink collision in the {q^^^ = 1, g'^' = 0, 
= 1) case. The solid line shows the initial energy density 
profile and the dotted line - the final. After the collision, 
there are four q ~ 2 kinks arranged so that the left most and 
the right most pairs of kinks form (g'^' — 2, g'^' = 2, g^"^-* — 
1) combinations and the two inner kinks form a (g'^-* = 2, 
g(2) _ 2^ — 0) combination. The rest of the original 
energy is being radiated away. The arrows show directions 
and approximate relative magnitudes of kink velocities. 



non-topological kink, which keeps propagating to the left 
unperturbed. Another way to describe this interaction 
is that the q = 2 wall has met its g = 2 "refiection" in 
the q — wall and formed a complex with it, while the 
remainder wall is radiated away. 

The (g(-^' — 1, g(^' — 1, g(^' = 0) case has essentially 
the same set of possible outcomes as the (q^^^ = 0, q^^^ — 
0, q'^^^ — 0) and will not be considered here. 

Fig. II illustrates the outcome of the kink-antikink col- 
lision with the initial (g*^^) = 1, q'^^^ = 1, g^^^ = 1) config- 
uration. The parameters were chosen to be h/X = —0.1 
and Vinitial — 0.2. The final configuration is that of two 
g = 2 kinks, arranged in a (g^^-* — 2, g^^' ~ 2, g^"^^ — 1) 
combination, moving away from each other. 

The outcome of a (g^^^ = 1, g^^^ = 1, g'^^ — 2) collision 
with h/X ~ —0.1 and Viniuai — 0.2 is shown in Fig. 0. 
The final configuration is that two (g*^^^ — 2, g^^^ — 2, 
— 1) combinations moving away from each other. 

The (g(i) = 1, g^^) = 2, g(3) = 1) case with h/X = -0.1 
and Vinitial =0.1 is illustrated in Fig. |[ Just before the 
collision, the g = 1 kink (originally on the left) collapses 
into three g = 2 kinks, with two outer kinks having large 
kinetic energies. Thus, an intermediate configuration is 
that of four g = 2 kinks, including the initial one (origi- 
nally on the right). The two right most g = 2 kinks are in 
a (g^^-* = 2, g(^) = 2, g'^' = 0) combination and, depend- 
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Figure 5: Kink-antikink collision in the {q^^^ = 0, g'^' = 2, 
_ case. Resulting configuration (dotted line) is that 
of a non-topological q = 1 kink moving to the left, while 
the q = 2 kink has captured its "mirror image" (the q = 2 
antikink) originally contained in the incident q = 1 wall. 



Figure 7: The {q'-^^ = 1, g'^' = 1, g'^' = 2) case with 
h/\ = —0.1 and Vi„itiai = 0.2. The final configuration is that 
two (g'^-* = 2, q^^'' = 2, g'^' = 1) combinations moving away 
from each other. 




Figure 6: The (g'^> = 1, g 



J2) 



= 1, 



J3) 



ft/A = —0.1 and 

of two q = 2 kinks, arranged in a (g 
combination, moving away from each other. 



= 1) case with 



0.2. The final configuration is that 

1) 



(1) 



^(3) 



ing on the initial velocity, may annihilate or chase each 
other forever. The two left most walls are in a (g^^^ = 2, 
— 2, q^"^) — 1) arrangement. 
A (7^^' — 1 kink and a g^^-' = 2 antikink arranged in 
a (q^^^ = 1, q^^^ = 2, (7^^^ = 2) combination repel. For 
h/X = —0.1 and Viniuai = 0.15 they scatter elastically. 



Figure 8: The (g^^^ = 1, g'^' = 2, g'^^ = 1) case with 
ft/A = —0.1 and Vi„itiai = 0.1. There are four g = 2 kinks left 
as a result of the collision. 



The lighter q = 2 kink (initially on the right) bounces off 
the heavier q — 1 kink and slows it down. 

The remaining two combinations from Table | are 

(g(i) = 2, ^ 2, ^ q) and (g(i) = 2, ^ 2, 



— 1). Since these are the only two initial com- 
binations that involve stable kinks and antikinks, they 
would be the most important ones if one studied evolu- 
tion of domain wall networks after the formation. How- 
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Figure 9: The elastic scattering in the (g'^' = 1, g'^-* = 2, 
q(^) = 2) case with h/X = —0.1 and ViniUai = 0.1. 



ever, these two are also the most "uninteresting" combi- 
nations from the novelty point of veiw. In the case of 
(q^^) = 2, = 2, g^^^ = 0), the dynamics and possible 
outcomes are qualitatively identical to the case of kink- 
antikink collision in the simple model. Namely, we 
observe a dependence on the incident velocity similar to 
that found in Refs. |, |, |, |, |l^, |l|| . The difference is 
that in the case of SU (5) the value of the incident veloc- 
ity leading to a given outcome depends on the parameter 
h/\ of the potential given in Eq. (|33|). 

In the = 2, q^^'^ = 2, q^^) = i) case, the kink and 
antikink repel and simply bounce off each other elasti- 
cally. 



VI. SUMMARY AND DISCUSSION 

As we have illustrated in the previous sections, kink- 
antikink interactions in SU{N) x Z2 can be of differ- 
ent types with many possible outcomes depending on the 
choice of the parameter values. 

In order to see how a pair of kinks will interact, one 
generally needs to specify the potential V{^) and evalu- 
ate the interaction energy between the two kinks. How- 
ever, we have shown that to a very good approximation 

the nature of the interaction can be determined by eval- 
uating Tr[Q(i)Q(2)]^ .^yjjgj.g Q(i) g^j^^ Q(2) 

are the charges 

of the two kinks defined by Eq. (p^. In particular, we 
have shown that kinks and antikinks may attract or re- 
pel depending on their relative orientation in the internal 
space. This is similar to interactions between global 0(3) 
monopoles, where the relative phase, not the topological 
charges, is what determines the nature of the interaction 



In our study of kink-antikink collisions in SU (5) x Z2 
we have seen a general tendency for larger mass topolog- 
ical kinks to split into fundamental kinks of the theory, 
such as the q — 2 kink. This suggests that kink-antikink 
interactions in SU{N) x Z2 can be described in terms of 
interactions between fundamental {i.e. globally stable) 
q = {N — l)/2 kinks. Namely, given the configuration of 
two kinks, one would look for the least energy configu- 
ration oi q = n kinks that would have the same global 
topology as the original pair of kinks. The outcome of 
the interaction would then be reduced to interactions be- 
tween attractive or repulsive pairs of q — n kinks. 

We have not investigated the detailed dependence of 
the dynamics of kink-antikink collisions on the initial ve- 
locities. Part of the reason is that the outcome strongly 
depends on the stability properties of each of the solitons 
and our numerical methods do not give us the possibility 
of properly accounting for all instabilities that can occur 
in the model. Each kink, except for the fundamental one, 
is unstable in a different way, namely, along a different 
direction in the internal space of SU(5), and also depends 
on the particular choice of the potential. It is possible 
that a more detailed study would reveal a connection 
with earlier work on Q balls |jl^ and global U{1) strings 
||l9| , where is was observed that at very high collision ve- 
locities fragmentation of solitons is generally suppressed. 
Such a study could be a subject of future work. 

From the equations of motion it follows that if the 
Higgs field components along non-diagonal generators of 
SU{5) were zero at the initial time, they would remain 
zero at all times, which was the case in our simulations. 
However, except for the q — 2 topological kinks, kink so- 
lutions can be unstable against perturbations along non- 
diagonal generators of SU{5) [|l3|. In a realistic domain 
wall formation scenario one would have to allow for all 
components of the Higgs to be excited and only stable 
walls would survive. Depending on the rate of the phase 
transition, temporary formation of unstable kinks will or 
will not be relevant. Nevertheless, stable and unstable 
kinks are solutions of the classical field equations and 
their interactions may become important whenever an 
exact or approximate SU{N) x Z2 symmetry is present 
in a theory. 
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